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Abstract 

We consider multidimensional gravitational models with a nonlinear scalar curvature term and form fields 
in the action functional. In our scenario it is assumed that the higher dimensional spacetime undergoes a 
spontaneous compactification to a warped product manifold. Particular attention is paid to models with 
quadratic scalar curvature terms and a Freund-Rubin-like ansatz for solitonic form fields. It is shown that 
for certain parameter ranges the extra dimensions are stabilized. In particular, stabilization is possible for 
any sign of the internal space curvature, the bulk cosmological constant and of the effective four-dimensional 
cosmological constant. Moreover, the effective cosmological constant can satisfy the observable limit on the 
dark energy density. Finally, we discuss the restrictions on the parameters of the considered nonlinear models 
and how they follow from the connection between the D-dimensional and the four-dimensional fundamental 
mass scales. 

PACS numbers: 04.50. +h, 11.25.Mj, 98.80.Jk 



1 Introduction 

Two of the most intriguing problems of modern cosmology are the problem of additional dimensions and the 
cosmological constant problem (CCP). The first problem follows from theories which unify different fundamental 
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interactions with gravity, such as M/string theory £Q, and which have their most consistent formulation in 
spacctimes with more than four dimensions. The problem can be naturally formulated as the following question: 
if we live in a multidimensional spacetime, why do we not observe the extra dimensions? Within the "old" 
Kaluza-Klein (KK) framework and the early Eg x i?s-heterotic string phenomenology the question is answered 
by assuming the extra dimensions so small (i.e. with a characteristic size r between the Planck and the Fermi 
scales 10 _33 cm < r < 10 -17 cm) that they are not accessable by present-day collider experiments. New concepts 
with the possibility for rich phenomenology opened up with the uncovering of D— branes by Polchinski 2 in 
1995. In " brane- world" scenarios of the Universe the usual 4— dimensional physics with its SU (3) x SU(2) x U (1) 
standard model (SM) fields is localized on a 3— dimensional space-like hypersurface (our world-brane) whereas 
the gravitational field propagates in the whole (bulk) spacetime. Depending on the concrete scenario there 
are different types of masking of the additional dimensions. Whereas in Arkani-Hamcd-Dvali-Dimopoulos 
(ADD) models |S] the extra dimensions are curled up to sizes smaller than 10 _2 cm, so that they are 
in agreement with present table-top Cavendish- type tests of gravity [S], they can be infinite in the Randall- 
Sundrum II (RS II) [7j and the Dvali-Gabadadze-Porrati (DGP) [5j model. In the latter models the appearing 
four-dimensionality of low-energy physics is achieved by inducing appropriate effective gravitational potentials 
on the world-brane. Beside their interesting phenomenology, brane-world models provide a possible resolution 
of the hierarchy problem. In ADD- type models this is due to the connection between the Planck scale Mpu^ 
and the fundamental scale M^ +D ^ of the 4— dimensional and the (4+-D')-dimensional spacetime, respectively: 

Mp m ~ Vd>M$£' DI) . (1.1) 

Vd> denotes the volume of the compactificd D' extra dimensions. It was realized in EI E] that localizing 
the SM fields on a 3— brane allows to lower M*(4 + p/) down to the electroweak scale Mew ~ ITeV without 
contradiction with present observations. Therefore, the compactification scale of the internal space can be of 
order 

r~Vn /D ' ~ l(#- 17 cm. (1.2) 

With Mew ~ ITeV, physically acceptable values correspond to D' > 3 jH] (for astrophysical and cosmological 
bounds see e.g. experimental bounds from table-top Cavendish- type experiments are given in 6 ), and for 
D' = 3 one arrives at a sub-millimeter compactification scale r ~ 10 -5 mm of the internal space. If we shift 
M„(4 + £ ) /') to 30TeV, as suggested in [5] (see also EE]), then the D' = 2 case satisfies all aforementioned bounds 
and leads to r ~ 1CU 3 — 10~ 2 mm. In order to not exclude this D' = 2 with its largest possible compactification 
scale r, we assume that the fundamental scale M*u+e>') can be of order 30TeV. Additionally, the geometry in 
the ADD approach is assumed to be factorizable as in a standard Kaluza-Klein model. I.e., the topology is 
the direct product of a non-warped external spacetime manifold and internal space manifolds with warp factors 
which depend on the external coordinates 1 . 

According to observations, the internal space should be static or nearly static at least from the time of 
primordial nucleosynthesis, (otherwise the fundamental physical constants would vary, see e.g. |12l I13j ). This 
means that at the present evolutionary stage of the Universe the compactification scale of the internal space 
should either be stabilized and trapped at the minimum of some effective potential, or it should be slowly 
varying (similar to the slowly varying cosmological constant in the quintessence scenario E3)- In both cases, 
small fluctuations over stabilized or slowly varying compactification scales (conformal scales/geometrical moduli) 
are possible. 

Stabilization of extra dimensions (moduli stabilization) in models with large extra dimensions (ADD-type 
models) has been considered in a number of papers (see e.g., Refs. El E3 EH EH EB 1213 l^) 2 - In the 
corresponding approaches, a product topology of the (4 + D') — dimensional bulk spacetime was constructed 
from Einstein spaces with scale (warp) factors depending only on the coordinates of the external 4— dimensional 
component. As a consequence, the conformal excitations had the form of massive scalar fields living in the ex- 
ternal spacetime. Within the framework of multidimensional cosmological models (MCM) such excitations were 
investigated in [221 1231 12*3) where they were called gravitational excitons. Later, since the ADD compactification 
approach these geometrical moduli excitations are known as radions |B1 116). 

Most of the aforementioned papers are devoted to the stabilization of large extra dimension in theories with 
a linear multidimensional gravitational action. String theory suggests that the usual linear Einstein-Hilbert 
action should be extended with higher order nonlinear curvature terms. In a previous paper |25j we considered 

1 The M-theory inspired RS-scenarios 171 1111 use a non- factorizable geometry with D' = 1. Here, the 4— dimensional spacetime 
is warped with a factor O which depends on the extra dimension and Eq. 11.11 is modified as: Mpiu\ ~ Q~ 1 M^w- in our paper 
we concentrate on the factorizable geometry of ADD-type models. 

2 In most of these papers, moduli stabilization was considered without regard to the energy-momentum localized on the brane 
so that the dynamics of the multidimensional universe was mainly denned by the energy-momentum of the bulk matter. A brane 
matter contribution was taken into account, e.g., in 1211 . 
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a simplified model with multidimensional Lagrangian of the form L — f(R), where f(R) is an arbitrary smooth 
function of the scalar curvature. Without connection to stabilization of the extra-dimensions, such models 
(4— dimensional as well as multi-dimensional ones) were considered e.g. in Refs. |26l 1271 l2~5| . There, it was 
shown that the nonlinear models are equivalent to models with linear gravitational action plus a minimally 
coupled scalar field with self- interaction potential. In |25| . we advanced this equivalence towards investigating 
the stabilization problem for extra dimensions. Particular attention was paid to models with quadratic scalar 
curvature terms. It was shown that for certain parameter ranges, the extra dimensions are stabilized if the 
internal spaces have negative constant curvature. In this case, the 4-dimensional effective cosmological constant 
A e f f as well as the bulk cosmological constant become negative. As a consequence, the homogeneous and 
isotropic external space is asymptotically AdS4. Because the considered nonlinear model is a pure geometrical 
one (only with a bare cosmological constant Ad as an exotic matter source included) the equivalent linear model 
contains only a minimally coupled scalar field as bulk matter. The null energy condition (NEC) T a f,N a N b > 
for this field reads T ab N a N b — (N a d a (j)) 2 > (with N a future directed null vector) and is satisfied only 
marginally when the internal spaces are completely stabilized and the scalar field is frozen out. Moreover, the 
weak energy condition (WEC) T a bW a W b > (with W a future directed time-like vector) is violated in this 
case because the energy density p of the scalar field is negative definite p < 0. As a result, the aforementioned 
parameters (the internal space scalar curvatures, A^ and A e ff) are negative in the case of stabilized internal 
spaces (see also [1710(11 121) ). 

However, a negative cosmological constant leads to a deceleration of the Universe instead to an accelerated 
expansion, as recent observational data indicate. According to these data our Universe is dominated by a dark 
energy component with negative pressure. For example, from observations of the clusters of galaxies follows that 
the energy density of the matter components which can clump in various structures is significantly undercritical. 
But, the position of the first acoustic peak in the angular power spectrum of the cosmic microwave background 
radiation (CMB) implies that the Universe is, on large scales, nearly flat. In other words, the energy density in 
the Universe is very close to the critical value. Thus, there must exist a homogeneously distributed exotic (dark) 
energy component j2H] . This observation is in agreement with the conclusion following from the Hubble diagram 
of type la supernovae (SN-Ia) at high redshifts, which also indicate that our Universe currently undergoes an 
accelerated expansion. Under the assumption of flatness, using the data of the CMB anisotropy measurements, 
high redshift SN-Ia observations and from local cluster abundances, the authors of Ref. (3U] found a constraint 
on the equation of state parameter ljq = P/p < —0.85 at 68% of confidence level. They concluded that 
this result is in perfect agreement with the ujq = — 1 cosmological constant case and gives no support to a 
quintessential field scenario with ujq > — 1. Results obtained in [31| also favor ojq rs — 1 at the present epoch. 

In Ref. we already indicated that the effective cosmological constant can be shifted from negative values 
to positive ones by including into the nonlinear model matter fields which satisfy the NEC. In the present paper, 
we demonstrate this effect explicitly by endowing the extra dimensions with real-valued solitonic form fields 
|32| . Such fields naturally arise as Ramond-Ramond (RR) form fields in type II string theory and M-theory. 
Within a generalized Frcund-Rubin settin g |33| th eir influence on the evolutionary dynamics of the Universe 
has been considered, e.g., in Refs. |34l 1351 l3ol \'tf\ and due to its simplicity we adopt this ansatz here for the 
stability analysis of our nonlinear model. From Eqs. (|2.7|l . (|2.8|) below, it can be easily seen that the real- valued 
form fields satisfy the NEC as well as the WEC. However, the strong energy condition (SEC) is violated in 
our model by the cosmological constant 3 . The presence of two types of fields in our equivalent linear model — 
the minimally coupled scalar field (which satisfies the NEC only marginally and which can violate the WEC) 
and the form fields (which satisfy both of these conditions) — leads to a rich and interesting picture of stable 
configurations with various sign combinations for the allowed cosmological constants as well as for the constant 
curvatures of the internal space. Beside stability regions with negative 4-dimensional effective cosmological 
constant h e ff < the parameter space contains also regions with A. e ff > which can ensure an accelerated 
expansion of the Universe. 

As mentioned at the very beginning of the Introduction, there still remains the problem about the incredible 
smallness of the cosmological constant [211 EH] • Moreover, it is completely unclear why its energy density is 
comparable with the energy density of matter just at the present time (the cosmic coincidence problem). Modern 
reviews on the cosmological constant problem can be found for example in |29U41ll4*2] . In our paper we show that 
for stabilized internal spaces a small positive cosmological constant of the external (our) space can arise from 
ADD- and KK-type multidimensional models. We demonstrate that the smallness of the effective cosmological 
constant can follow from a natural parameter choice of the considered nonlinear ansatz. Unfortunately, the 
extremely small value of the observed cosmological constant requires a very strong fine tuning of the parameters. 

The paper is structured as follows. The general setup of our model is given in section [3 In section [21 
we make the geometry of the spacetime manifold explicit — endowing the internal space with the structure 
of a warped product of n factor spaces (due to spontaneous compactification) . Furthermore, we specify the 

3 For a critical discussion of the different ECs we refer to 1381 . 
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generalized Freund- Rubin ansatz for the form fields and perform a dimensional reduction of the action functional 
to a 4-dimensional effective theory with (n + 1) self-interacting minimally coupled scalar fields (section |3J). The 
stabilization of the extra dimensions is then reduced to the condition that the obtained effective potential for 
these fields should have a minimum. In sectional we present a detailed analysis of this problem for a model 
with one internal space. The main results are summarized and discussed in the concluding section 

2 General setup 

We consider a D = (4 + D') - dimensional nonlinear gravitational theory with action 



ZK D JM L JM i=l a i- v ' 



(2.1) 



where f(R) is an arbitrary smooth function with mass dimension 0(m 2 ) (m has the unit of mass) of the scalar 
curvature R — R{g] constructed from the D-dimensional metric g ab (a, b = 1, ... , D). 

4 = 87r/M*+f D ,) (2.2) 

denotes the D-dimensional gravitational constant (subsequently, we assume that M*u+d') ~ Mew)- In action 

(12. 1|) . = Fm]ni...qi , t = 1, • • • j n is an antisymmetric tensor field of rank di (a di— form field strength) with 
indices from an index set sr{\ = {mi : max(m^) — min(mj) = di}, where rrii,ni, . . . , € sm. For simplicity, 

we suppose that the index sets suj, suj of tensors F^ with i ^ j contain no common elements as well 
as no indices corresponding to the coordinates of the -Do - dimensional external spacetime (usually Dq = 4). 
Additionally, we assume that for the sum of the ranks holds X)"=i d% = F> — Dq := D'. Obviously, this 
model can be generalized to tensor configurations F^\ F^> with intersecting (overlapping) index sets. In this 
case explicit field configuration can be obtained, e.g., when the indices satisfy special overlapping rules |36j . 
Such a generalization is beyond the scope of the present paper. Furthermore, we assume in our subsequent 
considerations that the index sets m,,n.j, . . . , 7^ do not contain the coordinates of the external spacetime 
Mo and, hence, the field strengths F^ can be associated with a magnetic (solitonic) p— brane system located 
in the extra dimensions as discussed, e.g., in Refs. |321 1351 131)] . 

The equation of motion for the gravitational sector of (|2.1|l reads 

f'Rab - \fg ab - V a V b /' + g ab Uf' = n 2 D T ab [F, g] , (2.3) 

where a, b = 1, . . . , D, f = df /dR, R ab — R a b\g], R = R\g\- V a and □ denote the covariant derivative and 
the Laplacian with respect to the metric g ab 

ifel _ -^abTJ v7_ _ a f f\=~\ —ab' 



□ = D{g] = 5 a "V a V h = -^=d a y/\g\ g ab db) ■ (2.4) 



Eq. (|2.3(l can be rewritten in the form 



f'Gab + \g ab (Rf - /) - V Q V fc/ ' + g ab Uf< = n 2 D T ab [F, g] , (2.5) 

where G ab — R ab — ^R~g ab , and its trace 

(D-l)U.f = jf-fR + K 2 D T[F,g] (2.6) 

can be considered as a connection between R and /. The energy momentum tensor (EMT) T ab [F,g~] is defined 
in the standard way as 

T ab [F,g] ^ — ^ — * L = ^T ah F^,g] , (2.7) 



i=l 
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where 

-(-- 

For the trace of this tensor we obtain 



PW a — — l--7j,F W p(i)m l n i ...q l , > T ? (i)n 4 ...g 4 

> f ,-/ I \ o™ m i n i ...q i - L ' u 'i ± ant .. .q t ± b 

n 

T[F,g]=J2T[F^,g 



with 



T 



2dj-D {i] F (t)m<n 4 ... g< 

2(di!) "*"*■••« 



The field strengths F^ satisfy the equations of motion 

i?W"i«i-9i. g . = o ( v /j^|F« m '"' =0, 

Vlsl V J ' q * 

and the Bianchi identities 



i = 1, . . . ,n . 



F r (i) ,=0, i=l,...,n. 



Following Refs. j^HlEIEHlj we perform a conformal transformation 



.9ab = & g, 



with 



ft=[f'(R)] 



(i b 



(2.8) 
(2.9) 
(2.10) 

(2.11) 
(2.12) 

(2.13) 
(2.14) 



and reduce the nonlinear gravitational theory to a linear one with additional scalar field. This transformation 
is well defined for f'(R) > (concerning the case /' < see footnote |HJ). The equivalence of the theories can be 
easily proven with the help of the auxiliary formulae 



and 



□ = rr 2 [□ + (£>- 2)g ab a.- 1 CL, a d b ] <^=> □ = n 2 [□ - (d - 2) g ab n- 1 n, a d b ] , 

R ab = R ab + §^(/T 2 V Q /'V b /' - (fr'VaVbf ^L_g abi fyi n f 

R = {f) 2/{2 - D) {R+^(fT 2 9 ab dafd b f' 2§-l(/')- 1 D/'} . 



Defining the scalar </> by the relation 



dR V D - 1 

and making use of (|2.15() - (|2.17|l . Eqs. I|2.5|) and l|2.6(l can be rewritten as 

G ab = K 2 D T ab [F, (f>, g] + T ab [(/>, g] 

and 

□ 



1 



y/(D-l)(D-2) 

The EMTs read 



D 



eV( D -'« B - ! ) — / - f'R + 



1 



y/(D-l)(D-2) 



i 2 D T[F,<j>,g] 



T ab = <l>,*4>,b ~ \9ab9 mn <t>, m (l>,n ~ ~<?a& e V<^H— > ^ ( R f - f) 



and 



T afc [F, 0, .g] = ^ eVc-'iMX [F« , g 

i=l 

n 2dj-D , 

T [F, (j), g] = e Vc°-D(d-2> 9 T F w , 5 



(2.15) 
(2.16) 

(2.17) 

(2.18) 

(2.19) 
(2.20) 

(2.21) 
(2.22) 

(2.23) 
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where T ab [F®,g], T [F« ,g] are given by replacing g — > g in Eqs. 
strengths F^ are now raised and lowered with the metric g. 
The equations of motion l|2.11[l for F^ transform to 



, (Trull. The indices of the field 



| 5 | e V(M(o-i) p'(*)Tn i n i ...g j 



Q. 



l,...,n, 



(2.24) 



whereas the Bianchi identities (|2.12J1 do not change. 

It can be easily checked that Eqs. (|2.19(l . Q2.20|l and (|2.24(l are the equations of motion for the action 



S : 




2k d 



R[g] - g ah <t>, a <t>,b - 2t/(</>) 



n 

«?>£ 

i=l 



1 



e V(D-ii(D-!) jrW' 



(2.25) 



where 



i?(0)e^ - / (R{4>))) 



B := 



D 



y/(D-l)(D-2) 



(2.26) 



and Eq. (|2.18(l is used to express R as a function of <fi : i? = R(4>). The scalar field (/> is the result and the 
carrier of the curvature nonlinearity of the original theory 4 (12. Correspondingly Eq. (|2.2()(l has a two- fold 
interpretation. It is the equation of motion for the field <j> and at the same time it can be considered as constraint 
equation following from the reduction of the nonlinear theory 1|2.1|1 to the linear one Ij2.25|l . Furthermore, we 
note that in the linear theory l|2.25[) the form fields are non-minimally coupled with the nonlinearity field <f>. 
(A minimal coupling occurs only for a model with n = 1, d\ = Dq, where according to l|2.10|l the trace of the 
form field EMT vanishes.) A comparison of the action functional with (|2.22l) shows that the last term in (|2.25|l 

coincides with the expression for the energy density — Tq [F, cj>, g] of the solitonic form field (due to F^J, = 

by the definition of F«). 

Let us consider what happens if, in some way, the scalar field <f> tends asymptotically to a constant: cp — > (f>o 
[precisely this situation should hold when the internal space undergoes a (freezing) stabilization]. From Eq. 
(I2.18|) we see that in this limit the nonlinearity in 12.11) disappears: f(R) « c\{R — Rq) + f(Ro) = c\R + C2, 
where C\ := f'(Ro) = exp(A<fio), Rq = R(cf>o), and — C2/(2ci) plays the role of a cosmological constant. In 
the case of homogeneous and isotropic spacetime manifolds, linear purely geometrical theories with constant 
A— term necessarily imply an (A)dS geometry so that the manifolds are Einstein spaces. In our model, the 
additional form fields destroy this asymptotical behavior. Instead, we obtain from Eqs. (|2.(i|) and l|2.5l) 



R 



D c 2 



1 



D - 2 ex ciD 



i=l 



2di - D 
2{dJ) 



(2.27) 



and 



Rab 



1 c 2 



1 



D-2cx c x D-2" u ^ 2(df!) 

i — 1 v ' 

where the form field product 



p(i) p{i)m i n i ...q i 
ratrii ...qi 



(i) Ui-.-qi 



(2.28) 



(2.29) 



is defined with respect to the metric g. For a form field, which asymptotically tends to a constant, the scalar 
curvature and the Ricci tensor also approach constant values. But whereas R and R are asymptotically connected 
by the relation [see Eqs. l(2~T7jl . $TH% and (5Z£fl ] 



R 



R, 



(2.30) 



the Ricci tensor R a i, will not be proportional to the metric ~g ab and, hence, the space will not be Einsteinian. 
This is in obvious contrast to a nonlinear model of purely geometrical type |25| where the stabilization will 
result in an asymptotical (A)dS £) spacetime. 

4 Thus, for brevity, we shall refer to the field <j> as nonlinearity scalar field. 
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In the rest of the paper we consider for simplicity a toy model 5 with quadratic curvature term: 



f(R) = R + aR 2 - 2A D , 
where the parameter a has dimension 0(m~ 2 ). For this model we obtain 

/' (R) = 1 + 2aR = e A4 > — • 77- _!_ f, 

and 



R= h [t 



U{<f>) 



1 



1 

4a 



(e 



A<j> 



1 



2A; 



(2.31) 



(2.32) 



(2.33) 



The condition 6 /' > implies 1 + 2aR > 0. In the limit a — > the nonlinearity is switched off and_the linear 
theory is recovered. Correspondingly, it holds /' — > 1 with implication ci = 1, 0o — * so that also R — > i? (in 
accordance with Eq. (|2.30l) ) and J7(0 — ► 0) — > A^. The corresponding region of weak nonlinearity is defined by 
the condition aR = e A<i " - 1 « 1. 



3 Spontaneous compact ificat ion and dimensional reduction 

The simple block-orthogonal structure of the field strength F shows that there is a preferable scheme for a 
spontaneous compactification of the multidimensional spacetime manifold: each of the form fields F^ % ' is nested 
in its own (^-dimensional space Mi. Thus, the D-dimensional spacetime M can be endowed with the structure 
of a warped product manifold 

M — > M = M x Mi x . . . x M„ (3.1) 

with metric 

n 

9 = g ab (X)dX a ® dX b = 5 (0) + £ g (l) . (3.2) 

i=i 

The coordinates on the (Dq = do + 1) - dimensional manifold Mq (usually interpreted as our (Do = 4)- 
dimensional Universe) are denoted by x and the corresponding metric by 

9 i0) =9$(x)dx»®dx» . (3.3) 

Let the internal factor manifolds M, be di-dimensional warped Einstein spaces with warp factors exp(/3 (x)) 
and metrics _ 

= e^^h&U {Vi)dyr ® d vT , (3-4) 



i.e., 

Rn 

and 



7 W 



= ^Imlm i rn = 1, . . . , (3.5) 



A'di = i?i - ferf^ , (3.6) 



where k = 0, ±1. The scale fj sets the characteristic size of Mi (modulo the warp factor exp(/3 )) and can 
be interpreted as an effective scale factor of the compact Einstein space Mi with metric j l and corresponding 

5 For considerations on higher order corrections to the gravity sector of M/string theory we refer to 1281 1431 . 

6 Obviously, the conformal transformation 12.131 . 12.141 becomes singular when f'(R) vanishes. The transformation itself can 
be extended from the /' > branch to the /' < branch with the help of an ansatz |27| Q = \f'(R)\ 1/( - D ~ 2) and a corresponding 
redefinition of the nonlinearity field <j>: e A ^ = |/'|. As result, one obtains an action functional for the /' < branch which differs 
from action 12.251 for /' > in its total sign and in the relative signs of the single terms as well as in the potential U(<j>). Most 
important, for a fixed sign of the Einstein-Hilbert term the kinetic term of the nonlinearity field has the correct sign, whereas the 
kinetic terms of additional matte r fields (in our case the form fie lds) have the wrong relative sign. This leads to a set of equations of 
motions which differ from 12.191 . I2.2UI . For details we refer to 1271 . Unfortunately, the sign switch occurs for all additional matter 
fields simultaneously and it is not controllable for some selected fields separately. Otherwise, it could have provided a natural 
mechanism for the generation of a phantom energy component with equation of state parameter u>q < — 1 (and a corresponding 
super-acceleration of the observable Universe) in the sense of 1441 . 
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volume 7 

V dl = j d d >y^y^\~rt i=l,...,n, (3.7) 



where has dimension 0(m di ). 

We note that the specific metric ansatz ()3.2JI - 1(3.41) for the warped product of Einstein spaces results in a 
scalar curvature R which depends only on x: R\g] = R(x). Correspondingly, the nonlinearity field 4> is also a 
function only of a;: <f> = 4>{x). 

The conformally transformed metric ((2.13(1 reads 

n 

g = tfg = (e A *f (D ~ 2) g := g® + £ e^W 7 « 0-8) 



with 



® = l***) m D ~ 2) 9<$, (3.9) 
P^ + ^Z*. (3.10) 
For the field strengths F"> we choose a generalized Freund- Rubin ansatz (see also |34l 1351 1551 137| 1: 

= V2 y/b®\ e mini ... q J {i) {x), FW^-« = (V2/^W|^ ^-^/W^) , (3.11) 

where = e 2/3 *7^^ and s mini ... qi is the Levi-Civita symbol. We use conventions where for Riemann spaces 
holds £ mi n i ...q i = E mini '" 9i and s mini ,,, qi e mini "' qi — di !. It can be easily seen that the ansatz l|3.11|l satisfies Eq. 

(12.24(1 (because and / depend only on x and the VTT^I factors cancel). The Bianchi identities l|2.12|l reduce 
to the equations 

afa*(x)/W(x)' 



with solutions 



dx^ 



/W(x) = A ) (3.13) 



(3.12) 



a , 



where en := e' 3 are the scale factors of the internal spaces Mi after conformal transformation 1(3.8(1 and /j = 
const . We choose the scale factors a, dimensionless so that the constants ff have dimension 0(m 4+D ) and 
ttp/j 2 ~ Q(m 2 ) . With 1(3. 13|) the energy density of the solitonic form field, and correspondingly the last term 
in action 1(2.25(1 . reads 

-Tg[F,<f>,g] = e^-Dt— ) * F^^F® = £ e»- J )' * := p(/3» , (3.14) 

Z— 1 2—1 « 

where for real form fields ff > 0. Again we see that for models with n — 1 and di — Dq this energy density 

decouples from the nonlinearity scalar field <f>; pf/3 1 , <f>) > pift 1 )- 

The fact that (f>, [3 Z and p depend only on x allows us to perform a dimensional reduction of action 1(2.25(1 . 
Without loss of generality we set the compactification scales of the internal spaces Mi i = 1, . . . , n at present 
time at (3 l — 0. This means that at present time the total volume of the internal spaces is completely defined 
by the characteristic scale factors r, (see 1(3.7(1 and footnote Q)): 

n p n 

= n / dSvOT ~ ]>* > ( 3 - 15 ) 

i=l„, j=l 



7 The volume is well defined for positive curvature spaces (k = +1). For compact negative and zero curvature spaces (A; = — 1, 0), 
i.e. compact hyperbolic spaces (CHSs) Mi = H di /Ti and tori Tj = R 3 /Tj , we interpret this volume as scaled volume of the 
corresponding fundamental domain ("elementary cell") Vd i ~ rf l X Vpn(i) (see, e.g., 1451 and references therein). Here H di , R d i 
are hyperbolic and flat universal covering spaces, and Ti, Tj — appropriate discrete groups of isometrics. Furthermore, we assume 
for the scale factors of the metrics ■y 1 ~ r^ 1 with 7* scaled in such a way that Vpu(i) ~ Thus, the volume Vd- is mainly 

defined by rj. In all three cases (fc = ±1, 0), the limit — > 00 results in an effective decompactification of the internal space with 
Vd- — > 00. In accordance with Eq. 13.61 . this means that the positive and negative constant curvature spaces flatten: R^. — * 0. 
Clearly, for compact Ricci-flat spaces holds Rd = by definition and without relation to the compactification scale of the torus. 



where D' = D — Dq = Y^j—i is the number of extra dimensions and Vd> has dimension 0(m D ). After 
dimensional reduction action l|2.25[l reads 



s = ^ ^mY[^{R 



Mo 



i = l 



7 (0) 



n 

£ R [<? W ] e- 2 "' - 2U(<j>) - 2 K 2 D p(J3, cj>)) 



(3.16) 



where Gy = diSij — didj (i.j = 1, 



1) is the midisuperspace metric |46l 14 7| and 



Vjy 



(3.17) 



denotes the Dq— dimensional (4-dimensional) gravitational constant. If we take the electroweak scale Mew and 
the Planck scale Mpi^ as fundamental ones for D— dimensional (see Eq. J2210 and 4-dimensional spacetimes 
(kq = Sit/Mp^) respectively, then we reproduce Eqs. I|l.l(l and (|1.2(l . 

Action H3.16|l is written in the Brans - Dickc frame. Conformal transformation to the Einstein frame |22II23| 



g(0) 



(3.18) 



yields 



S = 



^2 / d D "x^)\{R [g^] -Gifg^d^d^ - S^d^d^ - 2U eff (P^)} . (3.19) 



M 



The tensor components of the midisuperspace metric (target space metric on ) Gy (i,j — l,...,n) , its 
inverse metric G* J and the effective potential are correspondingly 



1 --■ S v 

= d l S lj H d,dj , G n = — 

3 J Dn - 2 J ' di 



1 



2-D 



and 



Ueff(P,<f>)= II' 



1 " 

- ^ Re~^ + U(c/>) + k 2 dP (P,, 



(3.20) 



(3.21) 



where U((f>) and p(/3, 0) are defined by Eqs. I|2.33|l and (|3.14jl . 

A stable compactification of the internal spaces is ensured when the scale factors of the internal spaces (3 l 
are frozen at one of the minima of the effective potential U e ff. The value of the effective potential at the 
minimum plays the role of the effective Dq— dimensional cosmological constant: Ueff\- n = A-eff- Assuming 
for the frozen scale factors at present time /3 1 = 0, we obtain the non-zero components of the asymptotic Ricci 
tensor (|2.28jl as 



d + ^Klffil 2 ^ ) g m ^ 



where 



1 c 2 1 



D - 2 ci ci D - 2 



(3.22) 
(3.23) 

(3.24) 



and Qq — (e j4 ^ >0 ) 1 ^ D 2 \ Thus, the asymptotic multidimensional spacetime is built up from Einstein-space 
blocks, but is itself a non-Einsteinian space due to the additional term in (|3.23() . 
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4 Stabilization of the internal space 



Without loss of generality 8 , we consider in the present section a model with only one d\ -dimensional internal 
space. The corresponding action (|3.19[1 reads 



S 



2«q 



d D °x 



\~9 (0) \{ 



R 



M 



where 



and 



g(0) 



<P 



gW^d^dvtp - g^d^d v <p - 2U eff (<p, 0)} , (4.1) 



'dip -2)^ 



U e ff{ip,<j>) = <D-2KD - 2 ) 



Do 



(4.2) 



(4.3) 



The potential U (4>) of the nonlinearity scalar field is given by Eq. i|2.33[l and the energy density l|3.14|l of the 
solitonic form field reads 

For brevity of the notation, we introduce 



(4.4) 



a := 2 



di(.D - 2) 
2d!- D 
C '~ v /(£>-l)(£>-2) : 

so that the effective potential reads 



6 := 2 
ft 



di 



_ (D-2)(D -2) 

h D Jl 



Ueff = e hv 



-R ie aip + U{(j)) + he c ^e adllfi 



(4.5) 



(4.6) 



From (|4.5|l we see that a real- valued form field /i implies a non- negative h = k 2 d ff > 0. For the rest of the 
paper, we continue to work with dimensionless scalar fields ip, <f> instead of passing to canonical ones (modulo 
8ir): (p = (pMpi(4), <f> = cj)Mpi^ and U e ff = Mpiu\Ueff- The restoration of the correct dimensionality is 
obvious. 

In order to ensure a stabilization and asymptotical freezing of the internal space Mi , the effective potential 
should have a minimum with respect to the scalar field tp 



d V U eff\extr = ' 



(4.7) 



so that for a minimum position at tpo = (which corresponds to a compactification scale (3 1 = at present 
time) it should hold 

^y^-Ri = bU(4>) + (adi + b)he c4 '. (4.8) 

This formula shows that the potential U e f / (ip, <j)) must also have a minimum with respect to <fr, because without 
stabilization of <f> the right hand side remains a dynamical function whereas the left hand side is a constant. 
This second extremum condition 

W//| extr = (4.9) 

yields 



[d^U + hce 



c<t>l 



0. 



Additionally, the eigenvalues of the mass matrix of the coupled (ip, < 
effective potential at the minimum position, 



J 



(4.10) 

-field system, i.e. the Hessian of the 

(4.11) 



8 The difference between a general model with n > 1 internal spaces and the particular one with n = 1 consists in an additional 
diagonalization of the geometrical moduli excitations. 
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should be positive definite 



1 



Tr{J) ± ^Tr 2 (J) - 4det(J) 



> 



According to the Sylvester criterion this is equivalent to the conditions 

Jn > 0, J 22 > 0, det(J) > 0. 



(4.12) 



(4.13) 



From H4.11fl we see that in the special case of 9^/7 e // 



= the Hessian is diagonal and the excitation 



modes of the fields decouple. The eigenvalues of J coincide in this case with the masses squared of the scale 
factor excitations (gravitational excitons \22\) m 2 = rr? v and the excitations of the nonlinearity field m| = mi. 

Let us now analyze the stability conditions l|4.8[l . (|4.1Q(I and H4.13|) explicitly. For this purpose we introduce 
the auxiliary notations 

<h ■= <t>\extr > X := e Ma > 0, q := 8aA D (4.14) 
and rewrite the potentials U, U e ff and the derivatives of U e ff at a possible minimum position (ipo = 0, 4>q) as 



= U\ extr 

Ueff\ extr 
d V U eff\ extr 



d <t> U eff\extr 



^W Ue ff\extr 
9 U U eff\ extr 
d <t><f> U eff\ extr 



(X-lf+q 



1 

8a 

-^Ri + U (X) + hX^ 



-Ri + bU (X) + (dm + b)hX^=^ = , 



J-X^t^ [(2A - B)X 2 - 2(A - B)X -(q+ 1)B] + hcX^A^- = q 
Rx - b 2 U (X) + [{d x a) 2 - b 2 ] hX^A^ , 



ia 
a 2 -b 2 



= cdiahX d-2 j 

= — X~t^ \(2A - B) 2 X 2 - 2(A - B) 2 X + (q+ 1)B 2 } + c 2 hX^^ 
8a 



(4.15) 
(4.16) 
(4.17) 
(4.18) 

(4.19) 
(4.20) 
(4.21) 



(The constants A, B are defined in Eqs. 12. 1811 and l|2.26[) . respectively.) We see that, for fixed dimensions Dq 
and d±, the two equations l|4.17|l . (|4.18|) describe a 3— dimensional algebraic variety V C M. in the 5— dimensional 



parameter (moduli) space 9 M. 



l + 3 (a, Ad, Ri, h, X). On the variety, inequalities l|4.13|l of the Sylvester 



criterion define then the parameter region T C V of stable compactifications. A natural strategy for extracting 
detailed information about the location of this stability region would consist in solving (|4.18() for X with 
subsequent back-substitution of the found roots into the inequalities l|4.13[) and the equation 1)4. 17p . In the 
following consideration we restrict our attention to the three simplest nontrivial cases which are easy to handle 
analytically. 



4.1 Zero effective cosmological constant: A e jf = 

By definition, we have A e // = Ueff\ t so that in the particular case A e ff = Eq. (|4.1(j|) yields the additional 
constraint 

Ue ff \ extr = -\Ri + U (X) + hX^ = . (4.22) 
Combining this constraint with 1)4. 17[) we obtain the relation 

adi-D 2d^ 

Rx = 2d 1 hX^^ - —?—U^(x) (4.23) 
d\ — 1 

2d 1 -D 

which can be used to eliminate the hX °- 2 term from (|4.18|l . As result we arrive at a simple quadratic 
equation in X with physically sensible solutions 



e Mo = X=l ~ — —2 » d i > 2 - (4.24) 

\ 1 + q , dx = 2 . 



^The compactification scale (modulus) n of the internal space Mi enters V C M via curvature scalar Ri (see Eq. 13.61 ). 



11 



With the help of Eqs. I|4.23ll . (|4.24|) and repeated use of a substitution-elimination technique, the potential Uq 
and the second derivatives (|4.19|) - 1)4. 21[) of the effective potential can be rewritten in the simpler form 

U (X) = ^f— -X-t^(X-I) (4.25) 



and 



Jll = 




J22 = 




J12 = 


d l<f, U eff 



i 2 diU (X), (4.26) 
j^X-t^ [EX + 4(A) - 1)] , (4.27) 



" l,a U (X), (4.28) 



' di - 1 

where E = (D - 4) 2 + 4(di - 2) > for di > 2. For the determinant of the Hessian (|4.11jl we get 

det ( J ) = I j ~ TTtFt MX)X-^ [(d x - 2)X + 1] . (4.29) 
a(cti — 1) — 1 

With the equations 14.23)1 - 1)4.29)1 at hand, we are well prepared to explicitly describe the location of the 
stability region T in the parameter space A4. Let us start with relation (|4.23|l . From the non-negativity 
conditions h > and e A ^° = X > we immediately conclude that for stable spaces Mi it should hold R\ > 
and Uo(X) > 0. Furthermore, we see from the latter condition and the Sylvester criterion J22 > 0, det(J) > 
[applied to (|4.27|l and l|4.29(l ] that for internal spaces of dimension d\ > 2 the parameter a is restricted to 
positive values 10 a > (the limiting case a — ► we discuss below). Finally, we note that Eq. 1)4.25)1 together 
with Uq(X) > and a > implies X > 1 and, hence, we find for d\ > 2 and d\ = 2 from the roots 14.24|) : 
(? = 8oArj > and also A^> > 0. 

Summarizing the obtained restrictions, we can describe the part O of the parameter space M. where the 
stability region T of the variety V is located: 

Tcvnece = (Q>o,Afl>o,]Si>o,/i>o,x>i)c^. (4.30) 

It remains to clarify what happens in the various limiting cases when the parameters reach the boundary <90. 

(L.l.l) q — > +0: According to (|4.24|) . this limit implies X — > 1, <f>Q — > 0. Because of q = 8aA^i we have to 
distinguish the two cases a — * and Ad — > 0. In these limits we obtain Uo(X) — > Au and /To(A^) — > 0, 
respectively. 

(L.1.2) a — > +0,A£> 7^ 0: The case a — > describes the transition to a linear model. Here we have U(<j>) — > 
Aj3, i?i — > [2di/(di — 1)]Ad and Ajj — > (di — In this limit, the mass of the tp-field excitations tends 
to infinity m\ — * to^ — > J22 — > 00 and the field itself becomes frozen at the position <^>o — * 0. The 
stabilization of the internal space occurs for R\,h,Au > with the gravexciton masses m\ — > m 2 — > 
J n = 4[(Z?q — 2)/(D — 2)]Aj3. This is in accordance with the results of Ref. [22], where a linear model 
with monopole terms was considered. 

(L.1.3) A D -> 0,a ^ 0: D ue to jQH and (L.l.l) this limit implies X ->• 1, <j> and !7 (I) -> so that 
according to l|4.28|l the excitation masses mi, TO2 decouple (J12 — > 0) and the gravexciton mass vanishes 
m\ — > m 2 — * 0. Hence, the limit A/) — * is connected with a destabilization of the internal space M\. 
The mass of the nonlinearity field excitations remains finite m\ — > m? ~ 1/a for a > 0. 

(L.1.4) h — > +0,a 7^ 0: From (|4.23|l - (|4.25|l we see that this limit of a vanishing form field is connected with 
R\,Uq{X),Ad — > +0. Thus the excitations of the nonlinearity field decouple from gravexcitons (J12 — > 0). 
Simultaneously, because of i?i — > +0 =>■ ri 00, the internal space Mi undergoes a decompactification 
and due to Uq(X) — > the effective potential C/ e // becomes flat in the (^—direction (Jn — * m v — > 0). 
This means that the internal space destabilizes, whereas t7 e / / remains well behaved with respect to cj). 
These results completely confirm the conclusions of paper |25j for a nonlinear gravitational model without 
form fields where a stabilization is only possible for A e // < 0. 

Finally, we note that for a model with d\ = Dq (and, hence, a vanishing trace of the form field EMT) the 
excitations of nonlinearity field 4> decouple from the gravexcitons: J12 = because of c = in 1)4.50 . (|4.28|l . 

10 Obviously, a negative a would yield a maximum of the effective potential U e ff instead of a minimum and our model would 
become unstable with respect to the conformal excitations of the internal space. The condition a > is also required in other R 2 
models 1481 to ensure tachyon-free configurations. 
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4.2 Traceless EMT of the form field: d x = D 



The easy handling of nonlinear models with a traceless form-field EMT as well as of models with a two- 
dimensional internal space Mi is connected with the structure of Eq. (|4.18|) . For X > 0, a ^ we obtain from 
Eqs. 14.5|l . (|4.18(l and the definitions of A and B: 



so 



[(D - 4)X 2 + AX - (q + 1)D] + (2dt - D)hX 2dl/(D - 2) = 0. 



(4.31) 



This algebraic equation reduces to a simple quadratic equation in X either when the last term vanishes due to 
2d\ — D = d\ — Dq = (the case of a traceless form field EMT) or when its degree l(d\) := 2d\/(D — 2) equals 
0, 1 or 2. For Dq = 4 we have l(d\ = 0) = 0, l(d\ = 2) = 1 and l(d\ — > oo) — > 2 so that as only sensible model 
remains Dq = 4, d\ = 2. It will be the subject of subsection 14.31 

For d\ = Dq = 4we find as physically sensible solution of Eq. H4.31(l 



X 



8q- 1 



We use this solution as well as the extremum condition (|4.17|l to rewrite Eqs. 14.15(1 - (|4.17|l . I(4.19[l 

the simpler form 



(4.32) 
|OT|l in 



Uq(X) 

Ri 



— x-^ 3 (x 

16a 

\uq{X) 



1), 



h , 



Ju 

J22 
hi 



-)2 



Ueff 



exir 



r) 2 



e ff I extr 
Ueff\ extr 



\Uq{X) + h 

[9h - U (X)} 



—X- 1 ' Z {2X 



= 



(4.33) 
(4.34) 
(4.35) 
(4.36) 

(4.37) 
(4.38) 



Obviously, there is no mixing of the excitations of the nonlinearity field cf> with gravexcitons ( J12 = 0) in this 
case: m^, = Ju, m 2 ^ — J22- Further, we read off from J22 > 0, X > that stable internal spaces are again 
only possible for a > and from (|4.32|) and X > that q is restricted to the half- line q > — 1. Additional 
information can be extracted by combining the condition Ju > with relations (|4.34l) . 14.35fl . what gives 

16h > i?i > 16U Q (X)/9 > 8A e// . (4.39) 

For the realistic case of a positive effective cosmological constant we find according to 14.3411 , l|4.39|l the conditions 

A e// >0: h > Rx/16 > U (X)/9 > h/3 > 0, (4.40) 

and, hence, from H4.32fl . (|4.33|) also the implication Uq(X) > X > 1 q > 0. We therefore conclude 

that such configurations are only possible for internal spaces with positive scalar curvature R\ > and positive 
bulk cosmological constant Au > 0. 

Let us briefly comment on some limiting cases. 

(L.2.1) h — * +0: In this case we recover the result of 25 that stable configurations are only possible for 
Ri, Ad, A e f f < [see the inequality chain l|4.39|l ]. 

(L.2.2) a — > +0, Ad ^ 0: For this transition to the linear model with freezing of the nonlinearity field at <fio — > 
and diverging excitation mass m| — > 00, the stability sector C M. can be read off from H4.39[) via 
substitution Uq[X) — > An- 

(L.2.3) Ad — * 0,q ^ 0: The limit is connected with q — > 0, X — > 1, </)q — > 0, Uq(X) — > and we have 
to distinguish two special cases. For a nonvanishing form field strength h ^ according to l|4.36|) . 
(|4.37|l inequalities Jn, J22 > hold so that both excitation masses remain finite. For vanishing field 



strength h — > +0 we obtain R\ — > 0, Ju — > 0, mf 
destabilization/decompactification with r\ — > 00. 



and the internal space M\ undergoes a 
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4.3 Two-dimensional internal spaces: d\ — 2 

According to Eqs. (|4.18|) . 14. 3111 . the extremum condition d<j>U e ff\ t = for models with two-dimensional 
internal space Mi and Do = 4 can be reduced to a quadratic equation and, hence, allows for an easy analytical 
handling of the models. Introducing the notation 

z := Aah (4.41) 

Eq. (|PT|) reads 

A 2 -2{z-l)X -3(g+l) = (4.42) 

and has solutions 

X 1>2 =z-l± ^{z-iy +%+T). (4.43) 

Furthermore, Eq. (|4.31(1 can be used to simplify the elements of the Hessian J. Setting Dq = 4 and d\ = 2 
everywhere in 14.5fl . (|4.19|) - 14.21[) and eliminating q with the help of (|4.42|) we obtain 11 



Ju = 




1 extr 


Jti = 


d l<f> U eff\ 


extr 


Jl2 = 




extr 



-^X- 1 ' 2 (5z + 1 - X) , (4.44) 
oa 

-Lx-V^X-z+l), (4.45) 
10a 

' X- 1/2 z (4.46) 



2V5a 



as well as 



det(J) = Xr— (X 2 - 6zX + 8z 2 - 4z - l) . (4.47) 

60a z X 

It is now easy to describe the part of the parameter space Ai where the stability region T C V is located. We 
start by substituting the solutions (|4.43|) into J 22 . Taking into account that X > 0, we get from the condition 
J 22 > for Xt, X 2 : 

± y/(z-l) 2 + 3(q+l)/a > . (4.48) 

Thus, the roots X\ and X 2 correspond to a > and a < 0, respectively. Because of z = Aah, h > this leads 
to negative values for X 2 so that this root is unphysical and stable configurations are restricted to X = Xi(z, q) 
and a > 0. The limiting case a — » +0 will be considered separately below. 

Furthermore, we see from the structure of Eqs. I|4.44|l - (|4.47H that the Sylvester criterion selects a region 
®(z,x) from the (z, X)— plane which can be interpreted as the projection of the stability region T on this plane. 
Explicitly we have 

Ju > =>■ X < bz + 1 , (4.49) 
J 22 >0 =>> X > z - 1 , (4.50) 

det(J)>0 =► H < l Z + V ^¥ 1 ^' ( 4 - 51 ) 
v ' \ X > 3z - y/z 2 + Az + 1 , v ; 

where the inequalities l|4.51|l are easily derived from (|4.47|l by calculating the critical values X c (z) for which 
det(J[z, X c (z)]) = 0. The intersection Q( z .x) of the sectors defined by the conditions (|4.49|) - (|4.51l) and X > 0, 
z > is shown on Fig. ^ 

In order to obtain information about the values of q = 8a (and A^) which allow for a stable internal space 
Mi it proves convenient to map the region Q/ Zj x) yieL quadratic equation 14.42fl or its solution X\(z, q) from the 
(z, A)— plane on an equivalent region @( z ,q) of the (z, q)— plane. For this purpose it is sufficient to transform the 
inequalities (|4.49l) - (|4.51|) and X > 0, z > for A and z into an equivalent inequality set for q and z. Let us 
demonstrate the mapping, e.g., for inequality ggg) . Substituting X = X x (z, q) = z - 1 + \](z - l) 2 + 3(q + 1) 
into the equation for the critical line A = X c (z) — 5z + 1 and solving for q we obtain as image of this line 
X c (z) a corresponding critical curve q c (z) = z(5z + 6) on the (z, q)— plane. (The same curve can be obtained by 
substitution of X c (z) into the quadratic equation 1)4.42(1 .1 With the help of two test points Pi — (zi, qi > q c (zi)), 



The curvature term in Ju of Eq. 14.191 cancels because of a = b = 1 for Do = 4, di = 2. 
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0.2 0.4 0.6 0.8 1 

z-4ah 



Figure 1: Projection ®( Zt x) of the stability region T c V C M. on the (z,X)— plane (shaded areas). The two 
lines A e ff = and R\ — (given in Eqs. JJjSSJl) separate the stable regions with: [X > X(z,A e ff — 0) : 
(A e// >0, R 1 > 0)], [X(z,A eff = 0) > X > X(z,Ri = 0) : (A e// < 0, R x > 0)] and [X < X(z, Ri = 0) : 
(Aeff < 0, i?i < 0)]. 



P2 = (z2,<l2 < Qc(z2)) above and below the critical curve q c (z), e.g. Pi = (1,26), P2 — (2,0), it is then easily 
seen that Xi(z, q) < 5z + 1 maps into q < z(5z + 6). Applying the same technique to l|4.50|l . (|4.51|) we obtain 

J u >0 q<z{5z + 6), (4.52) 

J 22 >0 =>■ g> -l-i(z-l) 2 , (4.53) 



detfJ)>0 ^ g < -1 + [4^ 2 + 10z + 1 + 2(2z + l)y g + 4z + 1] /3 , 

det^j>U ^ g > _l + [4z 2 + 10* + 1 - 2(2z + l)Vz 2 + 4z + lj /3 . 1 j 



Additionally we find from X > 



q > -1 - i(z - l) 2 for z>l, 
q > -1 for < z < 1 



(4.55) 



The resulting intersection region 6( z , g ) of Eqs. (|4.52l) - l|4.55|l is depicted in Fig. |21 

Let us now turn to the scalar curvature R\ and the four-dimensional effective cosmological constant A e f f = 
U e ff\ extr . The structure of Eqs. I|4.15ll . (|4.16|) and I4.17|l suggests to consider R\ and A e ff as functions of 
(z,X,q, a). Eliminating q from Eqs. (|4.16|) . I|4.17|l (with the help of 14.4211 1 we obtain 

i?i = -!-A- 1/2 (A + 4z- 1) , (4.56) 
6a 



A 



eff 



—X- 1/2 (X-2z-l). (4.57) 



12a 



The graphics of the functions 



R 1 (z,X) = Q => X\ {Rl=a) = l-4z 
A eff (z,X) = => X\ (Aeff=0) =l + 2z (4.58) 

are included in Fig. ^ For completeness, we map them also on the (z, q)— plane. Following the same scheme as 
above we obtain 

R 1 {z,q) = => q\ (Rl=0) = 2z(4z - 3) , 0<z<l/4 

Ae//=0 <?| (Ae//=0) = 2* (4.59) 
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0.2 0.4 0.6 0.8 1 

z-4ah 



Figure 2: Projection Q( z . q ) of the stability region T C V C M. on the (z,q)— plane (shaded ar- 
eas). The two lines A e // = and Ri = (given in Eqs. (|4.59|l ) separate the stable regions with: 
[q> q(z,A eff = 0) : (A e// > 0, i? x > 0)], [g(2, A e // = 0) > g > q{z,R x = 0) : (A e// < 0, i?i > 0)] and 
[q<q(z,R 1 =Q) : (A e// < 0, i?i < 0)]. 



(4.60) 



and the correspondences 

i?i > for l>l-4z, g> 2z(Az - 3) , 
A e //>0 for > 1 + 2z , g > 2z . 

From Fig. ^ an d Fig. we see that the nonlinear model with two-dimensional internal space Mi allows for 
stable configurations only in the cases 



A e //>0 for i?i>0, 
A e// <0 for sign (i?i) = ±1,0 



(4.61) 



It contains no stable configurations with an accelerated expansion of the Universe (A e /y > 0) for internal spaces 
of negative or vanishing scalar curvature R±. 
Finally, we comment on some limiting cases. 

(L.3.1) h — > +0, q 0: According to Figs. ^ 121 this limit corresponds to a vanishing form field z — > +0, 
and a stabilization is possible in the case of R\ < 0. Furthermore, for z — > we can approximate 
X = Xi(z, q) w (v — 1)(1 + z/v) with v := + 3q > 1 and the masses of the normal excitation modes 
of the coupled ip — cf>— field system follow from l|4.12ll . (|4.44|l - l|4.46|l as 



ml 



6a 
1 

10a 



(t,-l)-Va 



v+^z + 0(z 2 ) 
V -z + 0(z 2 ) 



2v 



(4.62) 



In the special case z = we completely reproduce our earlier results [35] on nonlinear stabilized models 
without form fields [d x = 2): m£ =m\ = -U (X), rn\ = m\ = X~ X I 2 (X + l)/(10a). 



(L.3.2) a — > +0, Ad, h ^ 0: For this transition to a linear model we have as in (L.l.l) U(X) - 
freezing of the nonlinearity field at (f>o — ► 0, X — ► 1. Using the approximation 

X = Xi(z, g) = 1 + z + (3q - 2z)/A - (3q - 2zf /64 + z 2 /4 + 0(a 3 ) 

we obtain the excitation masses as 

o 9 or 1 - 2' 



Ad as well as a 
(4.63) 



rnr. 



in 



to; 



= 3h - A 



D 



■f 0(a) -> 
0(a) > 



oo 



(4.64) 
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so that the freezing is clearly seen from the diverging mass of the nonlinearity field. Additionally, we find 
from gjjgj, P37jl 



Ri 

A eff 



A D + 3h-- [27(h + A D ) 2 - 8h 2 ] + 0{a 2 ) 



A D -h 



-a[{h + A D ) 2 -AA 2 D ]+0{a 2 ), 



(4.65) 
(4.66) 



what in the special case of a vanishing effective cosmological constant A e / / = reproduces the results of 
Ref. 22 for a linear model with Freund- Rubin form field: h = An = i?i/4, m 2 = 2h. 

(L.3.3) Ad -* 0, a,/i ^ 0: In this case we have q — > 0. A substitution of the approximation X = X\(z,q) = 

z - 1 + [(z - l) 2 + 3] 1/2 +3[(z- l) 2 + 3] ~ 1/2 q/2 - 0(q 2 ) into the Hessian shows that there is no special 
behavior of the excitation masses connected with this limit. 

(L.3.4) Ao,h — ► 0,a ^ 0: From 14.441) it follows in this limit Jn — > 0, so that beside a decoupling of the 



excitations the gravexciton mass vanishes m\ 



and the internal space M\ destabilizes. This is 



in full agreement with [21] where a stabilization for h = requires Ap < 0. 

(L.3.5) Ri — ► 0: In the limit R\ — > one observes a regular behavior similar to (L.3.3). For parameter points 
near the line X (z) := X( Rl=0 )(z) = 1 - 4z, < z < 1/4 we find from Eq. l|4*3H| A = X (z) + 

6aRiX^ 2 + 0{a 2 R 2 ) so that the Hessian yields excitation masses of the form m 2 2 {R\ ~ 0) = m\ 2 (Ri — 
0, z) + <7i^{z)aR\ + 0{a 2 R\) with some regular coefficients and 



1,2 



(Ri 



1/2 



9 - 5A T ^4(3 - 5A ) 2 + 5(A - 1) : 



(4.67) 



The masses m 2 2 (Ri = 0,z) have finite values except at the limiting points Xq(z — > 1/4) — > (or 
4>o — > — oo) and Xo(a — ► 0) — » 1 (or 0o — > 0) where both or one of the masses diverge. We see that, with 
exception of the limiting point 12 Xo(a ^ 0) — * 1, there occurs no destabilization of the internal space 
M\ for vanishing scalar curvatures R\. Due to the smooth behavior of the excitation masses under the 
transition R\ — > we can identify this limit with a stable decompactification r\ — > oo of an internal space 
Mi with fixed topology. Clearly, in our local approach a stable decompactificd space with r\ — > oo is 
indistinguishable from a stabilized internal space which is Ricci-flat from the very beginning. 



5 Conclusions and discussion 

In the present paper we investigated multidimensional gravitational models with a non-Einsteinian form of the 
action. In particular, we assumed that the action is an arbitrary smooth function of the scalar curvature f{R). 
Additionally, the D-dimensional spacetime was endowed with solitonic form fields of generalized block-orthogonal 
Frcund-Rubin type. This bulk matter ansatz leads to a naturally factorized geometry and a spontaneous 
compactification can be associated with it. For the considered models, we concentrated on the stabilization 
problem for the extra dimensions. As technique we used a reduction of the nonlinear gravitational model to a 
linear one with an additional self-interacting scalar field (nonlinearity scalar field). The factorized geometry as 
well as the generalized Freund- Rubin ansatz for the solitonic form field allowed for a dimensional reduction of the 
considered models and a transition to the Einstein frame. As result, we obtained an effective four-dimensional 
model with nonlinearity scalar field and additional minimally coupled scalar fields which describe conformal 
excitations of the scale factors of the internal space. 

A detailed stability analysis was carried out for the three most simplest configurations of a model with one 
internal factor space Mi and a quadratic curvature term: f(R) = R + aR 2 — 2 Ad, where Ad plays the role of 
a D— dimensional bare (bulk) cosmological constant. These three configurations are characterized respectively 
by: 1) a vanishing four-dimensional effective cosmological constant A e ff, 2) a traceless form-field EMT, or 3) 
a {d\ = 2)— dimensional internal factor space Mi. For all three configurations, a stabilization of the internal 
space is only possible in the case of a non-negative nonlinearity parameter a > and a bulk cosmological 
constant A^> restricted by the condition q = 8aAo > —1. The transition (Ad — > 0,h — > 0) is connected with 
a decompactification (i?i — ► 0, r\ — > oo) of the internal space M\. At the same time, it leads to a flattening of 

12 According to 14.671 . the limit Xo(a ^ 0) — ► 1 gives m\ — » m 2 — > 0, ra\ — > m? — * l/(5a). On the other hand, Xo(a ^ 0) — * 1 
implies z = Aah — » and according to 14.431 also Ajj — » 0. Thus, the results of (L.3.1) and (L.3.4) can be used to reproduce the 
same behavior of the excitation masses via 14.621 . 
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the effective potential in the direction of the scale factor excitations and, hence, to a destabilization of Mi (for 
a similar limiting behavior see also Ref. |25|L 

From the three configurations, the model with the two-dimensional internal space shows the richest features. 
It allows for stable configurations in the cases (A e // > 0, R\ > 0) and (A e ff < 0, any sign of as well as for 
Ricci-flat internal spaces (i?i = 0) . Interestingly, the various stable configurations belong to a connected region 
in the parameter space M and one can smoothly pass from one type of configuration to another one, including 
a transition to stable Ricci-flat internal spaces which can be described as "stable decompactifications" : R\ — > 0, 
t*i — > oo . As pointed out in the Introduction, such a rich picture became possible due to the presence of the 
real-valued form fields which satisfy the NEC and the WEC and which compete with the nonlinearity scalar 
field. The latter satisfies the NEC only marginally and can violate the WEC. 

Interestingly, for {d\ = 2)— dimensional internal spaces there exist parameter configurations with a, Ap>, h, R\ > 
that can provide positive values of the effective four-dimensional cosmological constant A e ff > (see e.g. Eqs. 
(|4.6U|) . 14.6111 1. Thus, an accelerated expansion of the Universe seems possible in accordance with observational 
data. Let us assume that the values of the bulk cosmological constant Ap> and the form field strength h are set 
at some characteristic scale Ad ~ h ~ M 2 . Then we find_for the parameters q ~ 8aM 2 , z ~ 4aM 2 and, hence, 
q ~ 2z. The latter corresponds toX-^l + 2z-^l + 8aM 2 (see (^33) an d comparison with (|4.59|) shows that 
such configurations should yield an almost vanishing effective cosmological constant A e // ~ 0. With the help 
of Eq. H4.56JI we can estimate the scalar curvature R\ of the internal space as 

z ,,o 1 z 4M 2 . , 

R l ~-X- 1/2 = 5.1 

a a Vl + 2z Vl + SaM 2 

On the other hand, its value is connected with the fundamental scales Af*(4+d 1 ), Mpi^ by the relations (|1.1|) . 
(H3 and iJSSJ: 

As mentioned in the discussion after Eq. (|2.33l) , the value of X can be used as a measure of the nonlinearity 
of the orig inal model: aR = e^ - 1 = X - 1. Wc see that weakly nonlinear configurations correspond to 
X m 1, whereas X ^> 1 leads to a strongly nonlinear regime. With the help of l|5.1|l and (|5.2() we express this 
dimensionless nonlinearity parameter X in terms of the different scales contained in our model: 

x-i + arfP-iaf * W^fr (5.3) 

\Ai*(4+di)/ \ M *{A+di)J 

From 1)5. 3|) we see that setting M <~ M*(4_)_ ( j 1 ) we obtain X ^ 1 for ADD-type TeV— scale models whereas X ~ 1 
can only be achieved for standard KK models with Mpiu) ~ ^I*(A+di)- Stably compactified internal spaces in 
ADD-type models can be obtained within a weakly nonlinear regime X ~ 1 if the bulk cosmological constant 
A/j and the form field strength h are related with the fundamental scales as 

A D ~fc~M*, M~lM, (4Wl) (ta) 2/ r (5 . 4) 

For M*(4 +( y -^1 — 30 TeV and di = 2 this implies M ~ 10 -4 — 10 _1 eV. It is interesting to note that this 
mass scale is of the same order as the lowest possible supersymmetry breaking scale m ~ Mg USY /Mpiu\ in the 
minimal supersymmetric extension of the standard model (MSSM) |3SJ with Msusy ~ 1 TeV. 

Above we demonstrated that the assumption A/j ~ h can result in a small effective cosmological constant 
A e ff. Let us now estimate the relation between A e ff and A^ , h in more detail and compare it with the observable 
value 13 of A e ff ~ 10~ 123 Ap-(4). For simplicity, we will restrict our consideration to a weak nonlinearity regime 
with X « 1, a > where the approximations l|4.65|l and (|4.(i()|) of (L.3.2) are valid. From (|4.(i6|) we see that 
to ensure a sufficiently small A e // the bulk cosmological constant A/j and the field strength h of the solitonic 
form field should be connected by 

h=(l + e)A D . (5.5) 
The value of the small e we will estimate now. With the help of relations H4.65[) and 15.4J1 we find 

10- 123 A PK4) ~ A e// « eA D (l-6aA D )/2 (5.6) 
A D - M 2 - R 1 



3 In our normalization conventions holds c = h = 1 and Ap 



~ al Pl(4) ~ ^p;(4)' 



18 



and, hence, 

10-- ~ e ( ^t*>) (5.7) 

so that 

e-10~ 65 (5.8) 

for d\ = 2, M„( 4+dl ) ~ 30TeV. According to l|5.6[l this value of e is not sensitive to changes of the nonlinearity 
parameter a in a weakly nonlinear curvature regime. Thus, we arrive at the conclusion that the ADD scenario 
in its simplest extended version can provide a simultaneous stabilization of the extra dimensions together with 
an adjustment of the effective cosmological constant to its observed value only in the case of a strong fine tuning. 
Although the solitonic form fields of our model are located in the compactified extra dimensions, the tuning of 
their effective energy density h to the bulk cosmological constant A^: h = (1 + e)Ar> is of a similar type as 
the four-form-tuning discussed in Weinberg's no-go theorem 0U| f° r a resolution of the cosmological constant 
problem (CCP). A shifting of the CCP to a parameter fine tuning is a rather general feature of models with 
compactified additional dimensions and form fields 14 . In a slightly reshaped form it also appears in the recently 
proposed brane- world model with two-dimensional "football" -shaped large extra dimensions (see also |52|) 
where the adjustment of the on-brane cosmological constant is shifted to an adjustment of the parameters of the 
off-brane "football". A possible resolution of the CCP for similar higher dimensional models with form fields 
following from an M-theory setup was presented in Ref. [HS|- Proposals for a resolution of the CCP within 
other scenarios comprise various anthropic approaches |Sni03], shifting of the CCP to a singularity problem |55j . 
possible graviton compositeness |56) . a holographic approach |57j as well as non-local modifications of gravity 
|58j . However, there is still no satisfactory and comprehensive solution of the CCP. The problem will probably 
remain challenging the scientific community until a final understanding of quantum gravity will be achieved. 
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